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Earlipr^i^uljs-ft^fcoherent  propagation  of  sound  in  correlatedrandom  distributions  of  two- 
pariunrtaf'particles  of  radius  a  (with  minimum  separation  bpla  small  compared  to  wavelength 
X=  2v/k )  are  generalized  topbuiin  the  refractive  and  absorptive  terms  and  the  corresponding  A  4 
bulk  parameters  to  order  (kaf.  The  present  development  includes  higher  order  terms  of  the  earlier 
multiple  scattering  by  monopoles  and  dipoles,  as  well  as  scattering  and  multipole-coupling  effects  town 

through  quadrupole  terms.  The  correlation  aspects  are  determined  by  the  statistical  mechanics  jL 
radial  distribution  function / (R )  for  impenetrable  particles  of  diameter  b.  The  new  terms  for  slab  * 

scatterers  and  spheres  involve  the  integral  oifR  (first  moment),  or  of  R  /In  R  for  cylinders.  The 
new  packing  factor  is  evaluated  exactly  for  slabs  as  a  simple  algebraic  function  of  the  volume 
fraction  tv,  and  it  is  shown  that  the  bulk  index  of  refraction  reduces  to  that  of  one  particle  in  the 
limiting  1.  Similar  results  are  obtained  for  spheres  in  terms  of  the  Percus-Yevick  approximation  | 

and  the  unrealizable  limit  je—  1 .  " 


I _  O  f-f 

l"®*  L U  co  tjg 

rs  s  ^ 


co 


PACS  numbers:  43.20.Ri,  43.20.Bi 


INTRODUCTION  0  e  0  ^  \ 

We  apply  general  results  for  coherent  propagation  in 
pair-correlated  random  distributions1  of  particles  with  mini¬ 
mum  separation  ( b )  of  centers  small  compared  to  wavelength 
(A  =  2  v/k )  to  obtain  additional  terms  in  k  of  the  bulk  param¬ 
eters  ( CJB )  and  index  of  refraction  (i/2  =  C/B)  considered 
before.2  Using  y  for  either  C,  B,  or  17 2,  and  the  form 
y  —  yr+  iya  +  iys>  for  small  kb,  we  obtained2  results  for  the 
refractive  ( yr )  and  absorptive  (ya )  terms  that  were  explicitly 
independent  of  k,  and  corresponding  results  for  the  scatter¬ 
ing  ( y ’,)  loss  term  to  lowest  order  in  k.  The  explicit  approxi¬ 
mations  for  yr  and  ya  for  spheres,  cylinders,  and  slabs 
(m  =  3,  2,  1,  respectively)  depended  only  on  the  particles’ 
radius  or  half-width  (a),  their  acoustic  parameters,  and  on 
their  average  number  (p)  per  unit  volume:  They  exhibited 
the  statistical  aspect  of  the  problem  only  in  the  volume  frac¬ 
tion  w=pv  with  v  =  v(a)  as  the  volume  of  one  particle.  The 
corresponding  scattering  terms  y,  were  additionally  depen¬ 
dent  on  (ka)m  and  on  the  low-frequency  limit  of  the  structure 
factor  W)  with  W  =  w(b  /2a)m  as  the  volume  fraction  of 
impenetrable  statistical  particles  with  diameter  b>2a,  i.e.,  in 
general,  each  particle  was  visualized  as  having  an  acoustical¬ 
ly  transparent  coating  of  thickness  (6  /2)  —  a.  The  present 
paper  applies  the  general  theory1  to  derive  the  leading  A- 
dependent  terms  of  y,  and  ya ;  these  depend  explicitly  on 
(ka\2  and  a/b  for  all  cases,  and  on  appropriate  correlation 
integrals^ W).  We  parallel  a  recent  development  for  the 
simpler  one-parameter  optical  case.3 

The  correlation  aspects  of  the  distribution  we  consider 
are  determined  by  the  statistical-mechanics  radial  distribu¬ 
tion  function4/)/? )  for  impenetrable  particles,  and  are  exhib¬ 
ited  explicitly  as  simple  integrals  over  all  R  of  the  total  corre¬ 
lation  function  F = /—  1.  The  integrals  for  spherical  and 
slab  particles  are  of  the  form  f  FR"dR  (moments  of  F),  but 
cylindrical  particlesalsoinvolvef  RF {In  R  )dR.  ThesecanaU 
be  evaluated  numerically  from  existing  statistical  mechanics 


forms  or  approximations4-15  for  /  We  obtained  explicit 
closed  form  approximations  before2-16  for  the  integrals  that 
arise  in  the  Tf  set,  and  also  used  the  required  ^integral  for 
spheres  in  a  related  development17  for  large  kb\  for  slabs,  we 
obtain  both  the  W  and  JV  integrals  from  our  earlier  Laplace 
transformation13  of  the  exact  Zemike-Prins  result12  for  / 
For  cylinders,  we  may  use  the  virial  expansion  for  /  to  con¬ 
sider  some  of  the  properties  of  jV. 

In  the  following,  for  brevity,  we  use,  for  example  (1:113) 
to  indicate  Eq.  (113)  of  Ref.  1,  as  well  as  essentially  the  same 
notation  as  before.1-3  We  generalize  the  earlier2  multiple 
scattering  monopole  and  dipole  approximations  for  C,  B, 
and  7 f  by  including  scattering  and  multipole  coupling  to 
higher  orders  in  k,  as  well  as  quadruples  (for  spheres  and 
cylinders)  to  lowest  order.  For  slabs  and  b  =  2a  (minimum 
separation  of  slab  centers  equal  to  slab  thickness),  the  explic¬ 
it  approximations  for  the  bulk  values  (y)  reduce  to  the  single 
particle  values  (/)  if  w— »1,  as  required  from  physical  consid¬ 
erations:  Y—*Y  because  the  particles  occupy  all  space.  The 
limit  in— <•  1  is  not  realizable  for  identical  spheres,  and  we  take 
wSwj  ss 0.63,  with  wd  as  the  densest  random  packing  intro¬ 
duced  earlier16  to  define  the  amorphous  solid.  However,  our 
explicit  approximations  for  y  for  spheres  also  reduce  to  y1  as 
w  increases  to  1;  we  regard  this  as  consistent  with  the  ap¬ 
proximations  involved  in  scaled  particle5  and  Percus-Ye- 
vick4,6  statistical  mechanics  theory,  and  with  the  closure  ap¬ 
proximation  used  in  the  multiple  scattering  theory.1  For 
cylinders,  we  take  wSwd  s0.84  as  before. 18,2  Were  an  anal¬ 
ogous  closed  form  available  for  the  jV  integral  for  this  case, 
we  would  expect  the  corresponding  approximations  for  y  to 
show  the  same  behavior  for  the  nonrealizable  limit  in— *1. 

The  present  application  of  the  general  theory1  to  larger 
kb  than  before2  plus  the  recent  applications17  to  large  kb 
provide  simple  forms  which  explicitly  display  the  functional 
dependence  on  all  key  parameters  for  many  practical  appli¬ 
cations.  Thus,  in  these  ranges  of  kb,  elaborate  machine  com- 
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putations  are  no  longer  required,  and  the  results  help  delin¬ 
eate  the  fundamental  physical  processes. 


I.  GENERAL  CONSIDERATIONS 

For  a  slab-region  distribution  and  a  normally  incident 
wave  <f>e  ~  *"*  (the  excess  pressure)  we  write 

<j>  =  eikz,  k  =  2v/A  =  2m)'//.t,  (1) 

with  Tje  as  the  index  of  the  embedding  medium.  The  corre¬ 
sponding  bulk  coherent  propagation  coefficient  is  given  by 

K  =  k7)b/r)e  =  kr),  T}2  =  C  /B,  (2) 

where  in  the  simplest  cases,  C  is  the  relative  compressibility 
and  B  ~ 1  is  the  relative  mass  density.  The  values  are  to  be 
expressed  in  terms  of  p  and  F  for  pair  correlated  particles 
specified  by  their  isolated  scattering  amplitudes  g( r,z).  The 
normalization  for  g  is  such  that  for  lossless  particles 

—  Re  g{x,z)  =  -Reg  =  ^|g(r,z)|2,  (3) 

with  as  the  mean  over  all  directions  of  observation  r.  The 
corresponding  known19  scattering  coefficients  an  are  nor¬ 
malized  by  the  form 


«  =  Sa'-:  an=°n(C'JB’-,x),  x=ka,  (4) 

where 

C'  =  ^,  C-  1=#,  B'- l=f, 

^  n  '  7 


with  C '  and  B '  as  the  particle’s  relative  acoustic  parameters, 
such  that  Im  C'>0  and  Im  fi '  <0  to  account  for  energy 
losses;  and  28  are  the  parametric  contrasts.  In  addition  to 
the  dependence  on  the  relative  parameters  and  on  x  (the  nor¬ 
malized  radius  or  half-width),  the  coefficients  depend  on  the 
dimensionality  (m)  of  the  problem. 

We  obtain  results  for  the  bulk  relative  values 
y  =  ( CyB.r/2 1  in  the  form 

r=Yi  +Yc  +iy,  =rt+x2rc  +  ixmr„  (6) 

where/’,  (with/  =  l,c,s)  depends  on  m.  The  forms  for/1,  and 
r„  corresponding  to  multiple  scattering  by  monopoles  and 
dipoles  to  lowest  order  in  x  for  the  real  and  imaginary  parts 
of  a0  and  a,,  were  discussed  before2  in  detail.  Now  we  obtain 
F‘ 

From  Rayleigh’s  results  for  spherical  dipoles,20  the  first 
approximation  for  sparse  uncorrelated  distributions  corre¬ 
sponds  to 


where  the  order  for  c  (and  subsequent  sets)  corresponds  to 
m—  1 , 2, 3.  In  Ref.  1 8  ( 1975),  for  lossless  particles  with  small 
parametric  contrast,  we  multiplied  Im  rjR  by  the  statistical 
mechanics  packing  factor  YT  to  obtain  the  appropriate  tj, 
for  the  correlated  case.  Appropriate  forms  for  C,  and  2?,  for 
spheres  and  slabs  were  given  by  Maxwell,2 1  and  for  cylinders 
and  spheres  by  Rayleigh.22  We  obtained  y,  and  y,  from 
(2:59),  which  followed  from  either  the  full  multiple  scatter¬ 
ing  procedure  or  an  equivalent  slab  procedure  (interface  ap¬ 
proximation).1'2  For  yc  the  two  procedures  lead  to  different 


results,  as  illustrated  for  circular  cylinders  by  ( 1 : 1 0 1  )ff.  Now 
we  apply  the  multiple  scattering  procedure  to  obtain  rc  ex¬ 
plicitly  for  an  unbounded  distribution.  See  discussion  after 
(1:102)  and  (23:37),  as  well  as  after  (23:50)  which  considers 
the  impedance;  discrepancies  may  arise  for  measurements 
on  plane-bounded  distributions  because  of  transition  region 
boundary  layers. 

Rewriting  (2:30)  in  the  original  form  (l:89)ff,  we  use 

_  JsLzJl  =  £  />„  =:P0  +  A  +  ^>2.  (8) 

c 

where  P2  does  not  arise  for  slabs.  The  coefficients  P„  are 
determined  by  the  linear  system  of  algebraic  equations  (ma¬ 
trix  equation), 

Pn=anV2n(l  +  'Z^n.V-n-yP.)>  (9) 

with  given  in  terms  of  the  correlation  integrals  by 

(2:21)ff.  All  cases  we  consider  are  covered  by 


P0  =  a0(1+P0jr0  +  ^  +  ^), 


Jr*) 

v3  r 


P2  =  a2r,*(l+-& 

\  rr 

Introducing 


.  +  +  : 

V 


V  } 


(10) 


D  D  95^  D  9^ 

* 2  ,  *  12  ,  2: 


(11) 


with  d„  —  dB—dc„  asin(l :95>— ( 1 : 100),  we  write  the  corre¬ 
sponding  bulk  parameters  as 

—  =  R2^-~~po 


m(m  -f 
Jj R^Pnd 


(12) 


~Pi- 


—  (PoV’2  ——(V2 

2) V  01  Hr' 


cP.) 


m(m  + 

~mlsp')  +  P2’  t>'2  =  cvb'-  <13> 


The  sets  of  d ’s  are  expressed  in  terms  of  Bessel  functions  of 
arguments  ka  and  Ka  in  the  forms  ( 1 :96)— ( 1 : 100);  we  use  the 
forms  in  terms  of /„  for  spheres,  in  terms  of  J„  for  cylinders, 
and  in  terms  =  cos  and  f  x  =  sin  for  slabs.  See  cited 
equations  for  full  details,  and  (1 : 101 H 1 : 107)  for  illustrations 
based  on  cylinders  and  for  comparisons  with  the  interface 
approximations.  As  shown  in  (1:99)  and  (1:100),  the  d  c  and 
dB  sets  are  proportional  to  #  and  28 ,  respectively;  thus 
C  =  1  if  C '  =  l.andP  =  1  if  B'  =  1,  as  required  by  the  theo¬ 
rems  (2:13)  and  (2:14)  discussed  earlier23  in  detail. 

Using  (11),  we  see  that  the  equality  ij2(B  —  1) 
—  (C  —  1)  =  1  —  i/2  is  satisfied  by  the  rigorous  series  forms 
in  ( 1 3),  (12),  and  (8).  To  the  accuracy  indicated  for  (6),  we  may 
replace  rj2  -f  cPx  by  1  —  cP0  within  the  large  parentheses  of 
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the  approximations;  thus,  the  approximations  in  (13),  (12), 
and  (8)  satisfy  the  equality  to  the  required  accuracy. 

The  isolated  scattering  coefficients19  for  the  slab, 
a„  =  a'/(  1  —  a' )  with  a'  imaginary  for  lossless  scatterers, 
are  given  by 

a’0=ix(<6  -x2  T0)/(  1  -  x2  D0)  +  0(x5), 
T0=[V{l+V'2)  +  2@v’2]/6, 

D0  =  (t,'2-1-2V)/2;  (14) 

a\  =  -ix(08  -x2Ty)/[\+08  -x2Z>,], 

Ty  =  [08(\  +  j}'2)  +  2V]/6,  Dy  =  (08  +  V)/ 2.  (15) 

For  the  cylinder,  we  use  a0  =  a’0/[l  —  Oq)  and  a„ 
=  a'/(l  —  a'n/2)  in  terms  of 

a’0  =  /(trx2/4)(#  -  x2  T’o)^  1  -  x2  D0)  +  O  (x6), 
T0=m  l+V'2)  +  3y’2]/8, 

D0  =  (t)'2  -1  +  2  VL  )/A=D  5  +  2  <4L  /4,  (16) 

where  L  =  ln(2/jcc')  with  c'  =  1.781...; 
a[=  -  i{irx2/4)(3  -  x2  Tx)/[  1  +  08/2  -  x2  Dx  ]  +  O (x6), 
r,  =  [^(l+-7'2)  +  2^]/8, 

Z>,  =  [3#  -  2^(1  +  2L )  +  rf'2  -  1)]/16 

=D\  -4  08L/\b-  (17) 

=  -  iir{x/2)*08/{2  +  08)  +  O{x%  (18) 

For  the  sphere,  a„  =  a'„ / [  1  -  a'„ /(2 n  +  1 ) ] , 

a'0  =  i(x3/3){V  -  x2  r0)/(  1  -  x2  Z)0)  +  0(x\ 
To=[(6(\+rt‘2)  +  2  087)'2n)/\O, 

2)0  =  (77'2-l+2^)/6;  (19) 

=  -  /(x3 /3)(^  -  x2  r,)/(l  +  ^/3  -  x2  Z>.) 

+  O  (x7), 

r,  =  [^(i  +i/,2)  +  2^]/io, 

Z),=  [3^-5^+2(V2-l)]/30;  (20) 

aj  =  -  ixs 2  08/9(5  +  2  08)  +  0(x7).  (21) 

Although  alternative  decompositions  may  be  con¬ 
structed  to  the  same  accuracy,  the  present  sets  of  T ’s  and  D ’s 
help  delineate  the  dependence  on  /§  =  C'  —  1  and 
08  =  B '  —  1;  comparison  of  a  given  set  for  m=  1, 2,  3  indi¬ 
cates  the  role  of  dimensionality.  If  08  =  0,  then  the  mono¬ 
pole  a0  dominates;  if  '€  =0,  then  the  dipole  a,  dominates. 
The  ratios  of  parametric  contrasts  in  (12)  and  (13)  do  not 
introduce  singularities;  the  explicit  approximations  dis¬ 
played  in  the  following  show  that  the  ratios  insure  fulfill¬ 
ment  of  theorems  (2:13)  and  (2: 14),  so  that  the  development  is 
consistent  to  the  required  accuracy.  For  m  =  2  and  3,  cylin¬ 
ders  and  spheres,  we  also  consider  rj2  for  rigid  particles 
(C'-*0,  B  '-*0)  as  the  limit  for  —  1  and  08-*  —  1. 

Using  (8)— (2 1 )  to  lowest  orders  in  x  of  the  real  and  imagi¬ 
nary  parts  of  all  components,  we  reconstruct  (2:63),  (2:67), 
and  (2:70)  in  the  present  forms  yx  and  Y,  corresponding  to  (6) 
and  use  the  results  to  simplify  the  subsequent  forms  for 
yc  —x2rc.  We  have 


C,  =  1  +  w # ; 


J9i  =  1  + 


w  08 
3  ’ 


m 


and 


y\  =  Cx/By  =  (1  +  w  #)/(!  +  w  08/3), 


(22) 


(23) 


where  /%  =  C‘  —  \  and  08  =  B '  —  1  are  complex  if  absorp¬ 
tion  is  present.  In  distinction  to  the  notation  in  Ref.  2,  the 
present  set  y,  includes  absorption  losses,  such  that 
Im  Cy  >  0,  Im  By  <  0,  and  Im^j  =  2  Im  tjx  Re  rjx  >  0.  The 
corresponding  scattering  terms  are 


C,  =x™ dmV2w7Tm,  Bs  =  — xm w  W ' ; 

m  32 

dm  ={l,?r/4,i>,  (24) 

with  packing  factors  as  in  (2:33), 


=■ 


(1  -  W) 


tn  +  1 


[1  +  (m-  l)wy 


(25) 


For  small  absorption,  i.e.,  for  |Re /|>|Im /|,  we  use 

#2cr|#|2  and  (08  /3)2?z\08  /3\2. 

From  (2:70) 


Cy  +  iC, 
By  +  iB, 


+  ivh 

082 
m3 2 


where  Re  (i;2/^,)  corresponds  to  attenuation  via  incoherent 
scattering  losses.  Similarly  for  the  present  generalization 
based  on  (6),  from  analogous  decomposition  of  rj2  =  C/B,  we 
obtain  the  additional  relation 


which  provides  a  check  for  the  independent  derivation  of  the 
new  terms  yc. 


II.  DISTRIBUTION  OF  SLABS 

For  slabs,  we  use  c  =  Hp/k  =  iw/x,  w  =  p2a,  and  re¬ 
tain  only  F0  and  F,  in  (8H 13)  with  =  c  +  and  ^ 

and  £0 ,  as  in  (1:177)  or  (2:23).  The  algebraic  system  for  P0 
and  Py  is  valid  for  all  ka  =  x,  but  we  consider  only  the  forms 
(1.4)  and  (15)  for  a„  and  the  corresponding  leading  terms  of 
the  correlation  integrals 

2?^2p  I*"  FdR  +  ik  2p[~  FRdR  -  1  -  ixN, 

Jo  Jo 

^y-bcijN,  (28) 

where  the  next  terms  are  O  (x2).  We  use  the  rigorous  func¬ 
tions3 

=  \\  —  Wf,  Ar=(l)»F(l-^+-5p), 

w=pb,  (29) 

with  W~  =  7Ty  of  (25).  For  b  =  2a  corresponding  to  W  =  w 
(minimum  separation  of  particle  centers  equal  to  particle 
width), 


31 


J.  Aooust.  Soc.  Am.,  Vol.  77,  No.  1 ,  January  1985 


Victor  Twersky:  Coherent  sound  in  scatterers 


31 


W  =  ( 1  -  w)2,  N  =  2w(l  —  Aw/3  +  W2/ 2).  (29') 

For  full  packing  u>— *•!  (the  limit  of  a  uniform  slab),  W-+0 
and  N-*  1/3.  Note  that  the  present  N  differs  in  sign  from  that 
used  earlier3;  all  packing  functions  are  now  defined  as  non¬ 
negative  for  w — ►!. 

Solving  (10)  to  the  required  accuracy  in  terms  of  the  7”s 
and/>’sof(14),  (15),  and  3)  =  1  -(-  SB(\  —w),  we  obtain 

cP0  =  —  wV%  —  x  2U0  +  ix^29SC); 

U0=T0-V  D0-  <g 2N - SB  Nrj2/®, 
T0-VD0=  [<g{l  -  i}'2  +  3V)  +  SB  t)'2]/3.  (30) 

CP,  =  rj2w{SB  -x2Ux  -  ix  SB27r/3))/S/>\ 

U,  -  [r,(l  +  SB)  -  SB  D,  +  SB2N}/3  +  /S  SB  N, 
T,(\  +  SB)-SB  D,  =  [SB(\- SB) +  <ig]/3.  (31) 

Substituting  into  (8)  gives  a  linear  equation  for  rj2  whose  so¬ 
lution  plus  (30)  and  (31)  also  determine  C  and  B  of  ( 12)  and 
(13).  To  lowest  orders  in  x,  we  obtain  the  sets  y,  =  7",  and 
y,  =  x r,  as  in  (22)-(26)  for  m  =  1 .  Thus  we  require  only  the 
k 2  correction  to  y„  i.e.,  yc  =  x2  rc  for  m  =  1. 

We  have 


2  x2w  (  Tr  ,  tV7i  ^  x2w  ( A  ,  A,  2\ 

3A0==‘#[i]'2  -  2  -  3#(1  -N)\- SB  t/'2, 

3A,=[SB(\-  SB  +3  SB  N)  +  <ig)/3)  +6<ig  SB  N, 

(32) 

9  iy'2C, 


:c=*2w(-t/0 


—  + 


SB  ) 


x2w<£ 


(y'2-2-3V(\-N) 


+  SB  rj'2w  +  rj\  (1  +  -Jf  N) 


)• 


-1L  + 

SB  j i2B, 

3)  + 

3 

+  33>) 

(33) 


x2wSB  (SB -\-3SBN+c<g[\-w) 

3  \  32 

+  ,■»«,),  (34, 

such  that  rfc  Bx  —  Cc  —  rj\  Bc  as  required  by  (27). 

If  W—  w — >1  (the  limit  of  a  uniform  slab),  then  lV-*l/3 
(as  well  as  £$-*\  and  y,— *Y\  and  the  set  yc  vanishes.  Since 
W-*Q  for  in-*’  1,  the  set  y,  also  vanishes,  and  the  results  for 
y  =  {rj2,  C,  B }  reduce  to  the  values  /  =  [if'2,  C',  B ')  for  a 
single  particle  as  required  by  elementary  physical  consider¬ 
ations.  * 

From  (23)  and  (26),  plus  (32)  with  rj\  replaced  by  rj2  x\ 
reconstruct  the  linear  equation  for  rj2  that  follows  directly 
from  (8)  in  the  form 

r?  «  1  +  w(V  +  x2A0  +  ix  V2?n 

~(w/3>)[SB  —  x2A,  —  ix[S92/3)7IS']ij2.  (35) 


The  first  approximation  is  j/2  =  as  in  (23),  and  iteration 
gives  Vi  =  W*  +  vt  +  ‘Vm-  We  proceed  similarly  in  the  fol¬ 
lowing  sections,  where  (8)  gives  a  quadratic  equation  for  ij2; 


in  order  to  minimize  elementary  algebraic  manipulations, 
the  yc  terms  are  listed  as  in  (32)— (34),  and  followed  by  the 
equation  for  rj 2  as  in  (35). 

The  results  simplify  for  the  special  cases  where  either 
one  of  the  relative  parameters  equals  unity.  Thus  if  B '  =  1, 
then  SB  =  0,  Bc  =  Bs  =  0,  and  B  reduces  to  unity  as  re¬ 
quired  by  the  general  theorem.23  For  this  case 

tj2  —  Cc  —  —  x2w  <i?2(2-w-  3N)/3,  (36) 

with  corresponding 

-q\  =  C,  =  l+w<if,  v]  =  C5=xw/r/S2.  (37) 


The  resulting  sum  C  =  C,  +  Cc  +  iCs  is  the  same  form  as 
(3:21)  for  e  in  terms  of  e'  and  8  =  e'  —  1. 

Similarly,  if  C '  =  1 ,  then  #  =  0,  and  C  reduces  to  uni¬ 
ty.  We  have 

Bc  =  -y2cB2t  =  x2w  SB2(2-w-3N)/332,  (38) 

with  corresponding 


B\  =  —■  =  1  + 

Vi 


Bs=  -V2B2 


SB  w 
3  ’ 

xw  W  SB2 

'  ~W2 


(39) 


For  the  present  one-dimensional  problem,  we  may  reduce 
the  forms  for  rj2  in  (38)  and  (39)  to  the  same  forms  as  in  (36) 
and  (37)  by  introducing  e'  =  l/B '  and  identifying  C  —  1  with 
# .  The  general  acoustic  two-parameter  results  (32)— (34)  also 
correspond  to  the  general  two-parameter  (e/i)  electromag¬ 
netic  cases  for  «|>  =  x<p;  if  <f*  =  E,  thenC'  =  t’  and  B '  =  1/fi'; 
if  <|*  =  H,  then  C '  =  p'  and  B '  =  1/e'. 

For  b  =  2a,  we  define 

=  w{2  —  w  —  3N)/2  =  w(2  ~7w  +  8 w2  —  3w*)/2 
=  w(2  —  3uj)(1  -  w)2/2,  (40) 


where  (1  —  w)2  =  W  decreases  monotone  from  1  to  0  as  u; 
increases  from  0  to  1.  As  discussed  before,  •¥ 
—  w  W  =  iu(l  —  w)2,  which  vanishes  at  u>  =  0  and  1,  has  a 
maximum  at  w—  1/3  corresponding  to  maximal  loss 
through  incoherent  scattering  (maximal  fluctuation  ef¬ 
fects).2,1®  The  function  SB /w  decreases  from  unity  at  w  =  0 
to  0  at  w  =  2/3,  reaches  its  minimum  value  at  w  =  7/9  and 
then  increases  to  0  at  w  =  1.  The  function  SB  which  vanishes 
at  w  =  0,  2/3,  and  1,  has  a  maximum  at  2w  =  1  —  (l/3)l/2 
and  a  minimum  at  2u>  =  1  -f  ( 1/3)1/2.  The  sign  of  the  correc¬ 
tion  terms  yc  changes  at  in  =  2/3. 


III.  DISTRIBUTION  OF  CYLINDERS 

For  cylinders,  we  use  c  =  iAp/k 2  =  iAw/irx2,  w  =  pwa2 
in  (8)— (13)  with 

2JTn  =c  +  <r0  +  ^2>  2^,2  =  cy 
2^22  =  c<172  +  1)  +  >n  +  ^4. 
in  terms  of  the  correlation  integrals  in  (1:71)  or  (2:22). 
We  have 

SBC 0~2irp  j  FRdR+  Up  j  Fln^y^R  dR 


=3T_  1+ijT, 


(FRdR= 

\  nir  /  J  nir 


(41) 
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where  the  next 


terms 


are  O  (x2 


We  work  with 


1  —  8  WFl ; 


w  vpb2 


F,  =  — 


du. 


(42) 


and 

-«=i-if.l»(A)+r-,„(^)  +  ,WT,: 


F,  =  -  f  F(ln  u)u  du.  (43) 

Jo 

Note  that  the  present  M  differs  in  sign  from  that  used  ear¬ 
lier.3 

We  may  evaluate  F,  and  F,  numerically  by  using  tabu¬ 
lated  values  of  F,  or  the  original  integral  equation  approxi¬ 
mations  in  the  computing  routine.9 

To  first  order  in  W,  we  use  the  virial  expansion14,3: 
F—  —  1  for  u  <  1, 


(44) 

and  F  =  0  for  u  >  2.  This  provides  results  for  W  and  M  cor¬ 
rect  to  0(lV2),  e.g., 


\  -  4W  +  y[3  \2W2/ir~  \  -  4W  +  6.6159W2. 

(45) 

The  closed  form  approximation  derived  earlier,18,2  i.e., 
5T=  2r2of(25), 

ar  =  (l  -  Wf/(\  +  W),  (46) 

gives  W  =  1  —  4 W  +  1W2  -f  ... .  For  the  unrealizable 
value  W=  1,  the  closed  form  W  vanishes;  as  shown  in  the 
following,  a  comparable  approximation  of  M  would  reduce 
to  3/4  for  b  =  2a  and  W=  1  in  order  for  y  to  equal  y1.  The 


corresponding  moments  are  then  F,  =  1/8  and 
F,  =  -  (In  2)/8  -  3/32. 

Solving  ( 10)  to  the  required  accuracy  in  terms  of  an  and 
the  F’sand  D’s  of  (16)— (18),  and  3)  =  1  -1-  3B{\  —  w)/2,  we 
obtain 


cP0  =  —  wi'iS  —  x2U0  +  iirx2  (€2  W/ 4); 
Uq  —  Tq—(€  D‘0 


+  ciS2M  /2  -  <€  3B(\  -  %r}r)2/43>, 
T0-VD^  =  [V(3-v'2)  +  3BV'2/*]. 
cPx  =  -  x2Ux  -  iir^SB2W/%9)/3)\ 


£/,  =  _  g  d  ; 

t  <g  39  [  1-7T) 


(47) 


jgy/l-y  JjKg, 

16  \  9  1 +B' 

\TX(2  +  3B)~  3B  D\  =  ^[3B(33B  +  4)4-4^].  (48) 
cP2  =  t}4wx2U2, 


u  =  &(i  +Bi)  _  £  ^  = 

2  8(1+F')  4(2  +  ^)l  +2^/  8^‘ 

(49) 

The  forms  for  t/2  illustrate  other  versions  that  also  arise  for 
Uv  Substituting  into  (8)  gives  a  quadratic  equation  for  rj 2 
whose  solution  plus  (47)-(49)  also  determine  C  and  B  of  (12) 
and  ( 1 3).  To  lowest  orders  in  x,  we  obtain  the  sets  y,  and  ys  as 
in  (22)-(27)  for  m  =  2,  and  as  before  we  require  only  the  set 
yc  =x2Tc. 

We  have 


V  l 


x2w  ( 

~bT\ 


-U0  + 


U>v  f 

9 


8/40  =  V(rj'2  -3-4  <6  M)- 9  v'2. 


x2w 


(Ao  +  ^-v2  t); 


164,  =  \3B{4  4-3  36  —4  3BM)  +  4<€\/3)  +  8  ^  <g(\  -  7T), 
1442  =  -^2(1  -  yr)/3)2-3B  [(1  +j9,)2/(1  +5')], 


(50) 


Cc  =x2«,(  -  C/0+  ^-g-C» ■■  + (7'2  -  3  -  4  3/+  0  77'2w  +  ^-[#(1  -  JT)  4-  1]),  (51) 


x?w3B  (3B  [ -3  +  4M  +  VHI- W\) -4  +  2V(l-w)  4V(1  -f1,,„l2n  ,  2(*,4-l)2,\ 

^ - 5 —  +  2’'i‘  +  ’,'lFTTi'> 


(52) 


The  set  (50)— (52)  satisfies  r\2c  Bi  =  Cc  —  ij\  Bc  as  required  by 
(27).  For  the  unrealizable  limit  W=m-*  1,  we  have 
—  3^*2— » 0  and  y,  =0,  y,  =  /;  if  we  use  M-+  3/4,  then 
yc-*0,  and  y— *y. 

The  original  equation  for  tj2  obtained  directly  from  (8) 
by  proceeding  as  for  (35)  is  a  quadratic  in  rj2, 


I 

v2  =  1  +  V2»^J 


JfL 


^  -x2A,  - 


iwx 


8  ^ 


4-  u»x2/42i74. 


(53) 
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The  first  approximation  gives  77? ,  and  we  obtain  1 j\  =  rj\ 
+  ifc  +  irj]  by  iteration. 

For  the  one-parameter  case  corresponding  to  5 '  =  1, 
we  have 

tj2c=Cc=  x 2w  &2(l+2w-4Af  )/8,  (54) 

ij]  =  C,  =  1  +  w  #,  tj2s  =  Cs  =  irx2u>  7T  #2/4. 

(55) 

The  resulting  sum  C  is  the  same  form  as  that  of  e  in  terms  of 
e1  given  in  (3:26)  for  the  optical  case  of  the  electric  field  paral¬ 
lel  (E|| )  to  the  axes  of  purely  dielectric  cylinders. 

Similarly  if  C'  =  1,  then 

Bc  =  -V\B2 

x2w  a  f3B'  + 1  -  a [4M+ 17?(1  -  arj] 

16  V  S2 

(56) 

B‘  "  (1+J )')/  1 

1+At, 

•.--M OT 

If  we  introduce  e’  =  1/5 '  and  <5  =  e'  —  1,  then  the  corre¬ 
sponding  sum  v2  =  V2i  +vl  +  ‘V2  =  €  is  the  same  as  in 
(3:33)  for  Ex.  The  general  acoustic  two-parameter  results 
(50)-(52)  correspond  to  the  two-parameter  (€41)  general  elec¬ 
tromagnetic  cases;  forE(( ,  we  take C '  =  e'and5 '  =  l/p‘;  for 
E1(  we  take  C'  —  p  and 5'  =  1/e'. 

For  comparison  with  (54)  and  (55)  to  0  we  have 

Bc~-v2c~- x2w  382(\  +  2w~4M+  2T)/ 16,  (58) 
5,-1  +  w  38  —  w{l  —  w)382/2, 

rj\~\ —  w  38 +w(\ +  w)382/2,  (59) 

5J  =  —  V2~  -  irx2 w  ar  ^  2/8.  (60) 

The  difference  in  the  statistical  factors  in  t])  of  (54)  and  (58)  is 
the  packing  term  but  both  factors  vanish  for  w  =  W-*  1 
and  M-+ 3/4. 

IV.  DISTRIBUTION  OF  SPHERES 

For  spheres,  we  usee  =  i4vp/k 3  =  iSw/x3,  w  =  pMra3/ 
3  in  (8)— (13)  with 

3^,,  =  c  +  JT0  +  2Jr2,  5 ^  =  207  +  2^  +  3^.,, 
35  JT22  =  c(7+  17i72)  +  7  JT0  +  10  MT2  +  18 

in  terms  of  the  correlation  integrals  ^  of  (1:148)  or  (2:2 1 ). 
We  have 

<ar0=  ar  _  1  +  /av*.  /x{2n  + 1),  (6i) 

where 

jr=l+4i rpj^  FR2dR  =  \  —  24fVF2; 
fF=^3 ”(y)*  ^  F(u)u2du,  (62) 


N  =  —  4npa 


FRdR 


-24(|K 


with  the  first  moment  as  in  (42).  The  next  terms  are  O  (x). 

By  differentiating  the  scaled  particle5  approximation 
for  the  equation  cf  state,  we  obtained 16  7T  =  of(25),  i.e., 

W  =  (\  -  W)4/(\  +2W)2,  (64) 

which  also  follows  from  the  Percus-Yevick  approximation.6 
The  corresponding  first  moment7  F,  of  the  Wertheim- 
Thiele  solution  of  the  PY  integral  equation6  gives 
e.vu  /  u/  u/2\ 


*=JfL_i”Wi—  E  +*L\ 

b  1  +2W  \  5  10  / 


For  b  =  2a  and  the  unrealizable  value  w  =  1,  it  follows  from 
the  closed  forms  (64)  and  (65)  that  W  =  0  and  N  =  9/5;  then 
as  shown  in  the  following,  y  —  Y-  See  (2:27)ff  for  virial  ex¬ 
pansions  in  powers  of  fV.  The  physically  realizable  range 
corresponds  to  W<.  Wd  =;0.63,  and  although  the  range  of  va¬ 
lidity  of  PY  theory  is  more  restricted,  the  limiting  behavior 
of  the  closed  forms  provides  a  useful  check.  The  simplicity  of 
(64)  and  (65)  obviates  detailed  computations,  and  their  rela¬ 
tions  to  the  rigorous  one-dimensional  results  (29)  emphasize 
the  essential  physics. 

Solving  ( 10)  to  the  required  accuracy  in  terms  of  the  an 
and  r,sand/),sof(19)-(21),  andS  =  1  +  38{\  —  u>)/3,  we 
obtain 

cP0  =  -w{V  -  x2V0  +  ix3#2  HT/ 3); 

0  0  0  3  9  S> 

T0-V  D0  =  h[V(4-V'2-5V)  +  38V'2],  (66) 

op,  =  77 2w[@  —  x2u  1  - ix3ds2  nrnsys)-, 

u  a 

V  3  9  /  9  9 


rjl3B22 

\JL  + 

2  w 

15 

3 

5 

r,[(3  +  ^)/3]  -  38  Dx  =  \[@(39  +  1)  +  #],  (67) 


2 

/3  +  25,' 

15 

1.3  +  25', 

We  proceed  as  before  for  (47)— (49).  We  have 

15^0  =  #  [77'2  -  4  -  5  <d(N-  1)J  -  38  rf'2, 

45 A,  =  {9 [38(38  +  1)  +  #]  -  5  382N}/&  +  10#  38  N, 

(69) 

3  3>2  5  (3  +  25')  1 


:,^v.+JL£Cl+l«) 

=  (x2u>  #/15)[t7'2  —  4  —  5  #(7V—  1)  +  , 
+  (77l/3^)(5^TV  +  9)], 
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D  2  (  Ux  ,  ^V'2Bt  ,  #  ,  2r,\ 

B'=xH-&+—t^+T¥+v'U2) 

x2w  3  (  -  3(^  +  !)  +  #(!  -w)  +  3  N( 5  +  2^)/3 

15  l 


5  ‘if  N 

3  3 


+  i7'25,  +  172 


2  2  (3  +  25 


5  (3  +  25 


0 


(71) 


The  set  (69)— (7 1)  satisfies  i)\  5,  =  Cc  —  if  Bc.  For  the  un¬ 
realizable  limit  W=  w—+ 1,  we  have  =  W2  =  0  and 
AT  =  9/5,  so  that  »/. 

The  original  equation  for  if-  obtained  from  (8)  by  pro¬ 
ceeding  as  for  (35)  is  a  quadratic  in  if, 

rj2=  1  +w[if  +  x2A0  +  (ix3/3)^2^] 


.J^3-x2Al-  5rV  + 

^>V  9  3  ) 

(72) 


corresponding  to  the  analog  of  (53)  for  m  =  3. 
For  the  one-parameter  case  5 '  =  1,  we  have 

if  =  Cc  =x2w  #2(6  +  3m  —  5iV)/15, 


(73) 


172  =  C,  =  1  +  m  tf,  ij2  =  Cs—  x2w  W  <if2/3.  (74) 
Similarly  if  C'  =  1,  then 

-35'  +  ^  Ar(5  +  2r/?)/3 


D  _  „2  DJ  _x2w^f 

- jj-J- 


32 


5,  ,  2  (3  +  25,)2\ 

d  f  ' 1  c  /->  .  ^  >v  y* 


5' 


5,= 


1 


1  + 


5  (3  +  25'). 
3  w 


5. 


Vi 

-V)B\  = 


3  ' 

-x3w  yr  32 


(75) 


(76) 


9  32 

where  3  =  1  +  .^(1  —  m)/3. 

If  we  introduce  1/5 '  =  V  and  1/5,  =  V„  then  the  re¬ 
sulting  form  of  if  equals 

•Tw 


tf  =  r,  =  i  + 


3. 


T  =V  -\, 


3 

3 


(77) 


Form  =  { 1,2,3}  we  have  3V  =  1  +  3^(1  —  m)(l  —  m  *)  as 
compared  to  3  =  1  +  3(1  —  w)m~'.  We  also  have  5,  3 
=  3JV,  3/Bx  3  =  -  r/^„,  so  that 

17?  -  -  5,/5  ?  =  *3m  2T  rV9  #*,  (78) 

bc  (-i  +  rN^  +  iify* 

{— 


31 


+  yl  +  Him+&\ 

11  5  3K'  +  2  / 


(79) 


For  the  simplest  acoustic  problems,  V  =  1/5 '  is  a  particle’s 
relative  mass  density.  The  electromagnetic  two-parameter 
and  one-parameter3  problems  for  spheres  give  different 
forms  than  (69H79). 

For  comparison  with  (73)  and  (74)  to  0(32),  we  have 
from  (75)  and  (76), 


Bez:-vl~  —  x2w  32(6  +  3m  —  57V  )7/(15)2,  (80) 

5,s:l  +  w3  -m(l  —  w)32/2, 
rj]  \  —  w  3  +  w(l  +  2u>)32/3,  (81) 

B, »  -  17? ss  -  x3w  W  32/ 9.  (82) 

Thus,  to  0(32)  the  statistical  packing  function  in  Bc  is  the 
same  as  in  Cc  of  (73)  and  (3:44). 

For  b  =  2a,  we  define 

^  _  m( 6  +  3m  —  5,#)  _  w( 2  —  5w  +  4 w2  —  w3) 

~  6  2(1  +  2  m) 

m(2-m)(l-m)2 
2(1+ 2u>) 

as  the  analog  of  (40).  The  relation  of  3  to  of  (25)  is 
indicated  by  (1  —  w)2/(  1  +  2 w)  —  Wxf2.  The  functions  W 
and  3 /w  decrease  monotone  from  1  to  0  as  m  increases  from 
0  to  1.  The  product2,18  y  —  wW~  has  a  maximum  at 
ms; 0.13,  and  3  has  a  maximum  at  ms; 0.22. 

V.  COLLECTIVE  FORMS 

The  results  (35),  (53),  and  (72)  correspond  to  (8),  i.e.,  to 
if  =  \  —c  2 Pm  expressed  as 

172  =  1  +  m(#  +  x2A0  +  ixmdm  <g 2  JT) 

(84) 

with  A2  =  0  for  m  —  1.  The  solution  of  the  matrix  equation 
(10)  for  the/*,  isolated  corresponding x2Un(i]2)/c  terms;  the 
regrouping  of  the  sum  of  U'%  in  (8)  in  order  to  make  the 
dependence  on  if  explicit,  introduced  the  coefficients  An. 
The  explicit  versions  are  special  cases  of 


+ 


rv 


5,  \  3 

ixmd„wW  > 


7?) 


+  ■ 


(85) 


5,  V  '  m32 
with  coefficients  defined  in  (32),  (50),  and  (69). 

Similarly  (33),  (51),  and  (70)  are  covered  by 

V  ,23Ci  Vy\  \ 

m(m  +  2)  (m  +  2)3  ) 

+  ixmdm  #2m  W,  (86) 

and  (34),  (52),  and  (71)  by 

B  =  B,  +  -  -g-  +  ^  U,) 

\  3  m(m  +  2)  (m  +  2)3  2) 


C  ssC,  +  x2w^  -U0  +  -2-— — -i  + . 


dm  3f 

m32' 


(87) 


with  the  Un  defined  in  (30),  (31),  (47H49),  and  (66H68),  for 
m  =  1,2,3. 

VI.  RIGID  PARTICLES 

We  derive  the  analogous  results  for  if  for  distributions 
of  rigid  spheres  and  cylinders  as  the  limits  of  the  two- param¬ 
eter  functions  for  C  '—>0  and  B  '—>0.  We  consider  spheres 
first,  because  explicit  closed  forms  in  IF  are  available  for  the 
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required  packing  functions  W and  N.  Then  we  consider  cyl¬ 
inders  and  use  the  development  to  infer  an  additional  prop¬ 
erty  of  the  implicit  function  M{  W). 

For  rigid  spheres  we  substitute  #  =  —  1,  =  —  1, 

and  3>  =  1-1-  SS{\  —  w)/ 3  =  (2  +  w)/ 3  in  (69)  to  obtain 

9-5  N  __  27-5A(l  +2u>) 

0  15  ’  '  45(2  +  w)  ’ 

21(10  -^-4SAf| 

[15(2  +  10)1’ 

Although  these  vanish  for  W  =  w  =  1  (for  which  case 
A  =  9/5),  we  show  that  for  W  —  w  each  coefficient  contains 
1  —  w  as  a  factor;  this  enables  us  to  use  the  same  iteration 
procedure  in  (72)  as  for  the  two-parameter  case  even  for  the 
unrealizable  limit  of  w—*l. 

Substituting  (88)  into  (72)  for  #  =  38  =  —  1,  we  have 


rj2  =  1  —  w  + 


hv 


+  ■ 


2  +  w 
iwx 3  W 


I)2  -|-  x2w\ 


1+- 


(ao 

)■ 


+ 


3  A, 

2  +  w 


V2  +  A: 


(89) 


j  \  (2  4-  w)2 

Solving  by  iteration  yields  the  leading  term  of  the  real  part 

rj2  =  (2  +  w)/2  =  1  +  w/2,  (90) 

as  well  as  the  leading  term  of  the  imaginary  part 
X3w  W[1  +  2 w) .  ^  _  (1  -  W)4 


_ 
V s  = 


12(1  —  w) 


(1+2W)2 


W>w, 

(91) 


(92) 


and  the  x2  correction  to  the  real  part 

2  w[2  +  w)  /  3_  (2  +  M2)2\ 

Vc  2(1  -w)  V  0  A'  2  +  4  / 

The  leading  terms  follow  directly  from  (23)  and  (26),  or  from 
the  closed  form  (2:74);  the  factor  1/3  in  (2:75)  should  be  re¬ 
placed  by  1/12.  The  physically  realizable  domain  corre¬ 
sponds  to  W<  Wd  =0.63,  and  w<  W.  For  W—  w,  (91)  has  a 
maximum  at  u>=0.16. 

It  is  clear  for  (91),  that  even  for  W  =  w— *1,  the  iteration 
procedure  is  valid:  1  —  w  is  a  factor  of  and  W /{ 1  —  u;) 
vanishes  as  0(1  —  in)3.  Similarly  for  (92),  for  W  =  w  in  N  of 
(65)  we  may  factor  1  —  w  from  the  coefficients: 

A„  =(1  -  w)s/n- 


S/  Q  - 


sf  5  — 


=  - 


3  —  w  +  w2 
5(1  +  2w\ 

2(20  +  2w  +  5  w2) 
45(2  +  uj)2(1  +  2w) 
Substituting  into  (92)  we  obtain 
2  _  x2w{2  4-  w ) 

Ve  ~  180 


9  —  w  +  w2 
15(2  +  w)  ’ 


(93) 


/74  -  16w  +  23  w2 

9(9  -  w  +  uM 

\  1  +  2  w 

2  +  w  ) 

jAnfl  —  w) 
180(1  +2u>) 


(67  —  44iu  —  5in2), 


(94) 


which  has  a  maximum  at  u>=0.28  and  vanishes  in  the  unrea¬ 
lizable  limit  We  use  these  limiting  results  as  guides  for 
the  analogous  case  of  cylinders. 

For  cylinders,  we  substitute  9%  =  £%  =■  —  1  and 
3>  —  1  +  38 \\  —  w)/ 2  =  (1  +  w)/2  in  (50)  to  obtain 


A0 
A2  = 


3-4 M 
8  ’ 
W 


2  (3 -43/) -8  ,  (1 

'4l_  8(1  +w)  +  2  ' 


(95) 


4(1  +  w)2’ 

which  vanish  for  W  =  w  =  1,  if  we  use  CM/'  =  W2  =  0  and 
M  =  3/4.  Substituting  into  (53)  yields 


2  ,  ,  2WJ]2  2 

t]  =  1  —  w  -| - —  +  x  w 

1  +  w 


U  +  ^t  +  A,v‘) 

\  1  +  w  ) 


toW  _2gM 
4  l  (1  +  w)2)’ 


(96) 


and  iteration  gives 
rj\  —  1  +  w, 

2  ttx2w  2T(3  +  w)  (1  -wf 

Vs  4(1  -w)  ’  l  +  W’ 

as  obtained  directly  from  (23)  and  (26),  or  from  the  closed 
form  (2:74);  see  (2:76).  In  addition,  the  correction  to  is 
given  by 


(97) 

(98) 


t)2c  =x2w  [(1  +  w)/{  1  —  m)]  Mo  +  2/4 1  +/42(1  +  ui)2]. 

(99) 

The  physically  realizable  domain  corresponds  to  W <Wd 
=0.84,  and  tn<  W. 

For  (98)  and  JV  =  in,  we  factor  1  —  w  from  7?~(w),  and 
the  result  for  rj]  vanishes  for  the  unrealizable  in  =  1 .  For  (99) 
however,  although  the  individual  An  vanish  for  the  limiting 
values  7T  =  0  and  M  —  3/4,  this  does  not  insure  that  rj2c 
vanishes.  To  eliminate  a  nonvanishing  residue,  we  require 
that  for  large  in, 


3  w  (1- m)8  +  C>(l -in)" 
3  +  in 


(100) 


with  n  >  1.  This  inference,  based  on  the  behavior  of  rjl  for 
spheres  in  terms  of  the  known  closed  forms  for  N  and  W 3, 
may  facilitate  development  of  a  corresponding  closed  form 
for  M  to  use  with 
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